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Impacts of a power-law non-thermal electron tail on the 
ionization and recombination rates 
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Abstract. We have investigated the effects of a non-thermal electron population on the ionization and recombina- 
tion rates. The considered electron distribution is defined as a Maxwellian function below a break energy £?b and 
as a power-law function of index a above this energy. We have calculated the collisional (direct and excitation 
autoionization) ionization coefficient rates as well as the (radiative and dielectronic) recombination rates. Practical 
methods are given to calculate these rates in order to be easily included in a computer code. The ionization rates 
are very sensitive to the non-thermal electron population and can be increased by several orders of magnitude 
depending on the temperature and parameters of the power-law function (Eb and a). The non-thermal electrons 
have a much weaker effect on the (radiative and dielectronic) recombination rates. We have determined the mean 
electric charge of elements C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni for different values of the break energy and 
power-law index. The ionization balance is affected significantly, whereas the effect is smaller in ionizing plasmas. 
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1. Introduction 

The ionization and recombination rates for astrophysical 
plasmas have usually been calculated for a Maxwellian 
electron distribution 



(e.g., Arnaud & Rothenflug 1985, 
Arnaud & Raymond 1992, Mazzotta et al. 1998] ). However, 
in many low-density astrophysical plasmas, electron distri- 
butions may differ from the Maxwellian distribution. The 
degree of ionization of a plasma depends on the shape of 
the electron distribution, as well as on the electronic tem- 
perature. This has been studied for the solar corona (e.g. 
Roussel-Dupre 1980[Owocki & Scudder 1983, Dzifcakova 



1992 , Dzifcakova 1998) a nd for evaporating interstellar 



clouds (Ballet et al. 1989 ), where a non-thermal electron 
distribution occurs in places where there are high gradi- 
ents of density or temperature. 

A non-thermal electron population is expected in var- 
ious astrophysical plasmas. Strong shocks can convert a 
large fraction of their energy into the acceleration of rel- 
ativistic particle s by t he diffusive shock ac celera tion pro- 
cess (e. g., Dr ury |l983[ Blan dford & Eichler |l987| , Jones & 
Ellison 1991 , Kang & Jones 1991 ). Direct evidence for the 



presence of accelerated electrons up to relativistic ener- 
gies (~ lGeV) comes from the observations of radio syn- 
chrotron emission in supernova remnants and in clusters of 
galaxies. More recently, non-thermal X-ray emission has 
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been reported in several shell- like supernova remnants and 
interpreted as synchrotron radiation fro m cosm ic-ray elec- 
100 T eV (K oyama et al. 



1995 



Allen et al. 



trons up to 

1997L Koyama et al. [19971 , Slane et al. |1999| , Slane et al. 

2rjoT|) . 

A number of recent works have focused on the non- 
thermal emission from supernova remnant s (e.g . , Laming 
20 01, Elli son et al. [2000[ Ber ezhko & Volk |2000j , Bykov et 
al. |2000b[ B aring et al. |1999| Gais ser et al. 1998j , Reynolds 



1996| , [1998| , Sturn er et al. |1997| ) and cluste rs of galax- 
ies (e.g., Sara zin 1999 , Bykov et al. 2000a , Sarazin & 
Kempner 199E ) . The impact of efficient acceleration on the 
hydrodynamics and thermal X-ray emission has been in- 
vestigated (Decourchelle, Ellison, & Ballet 2000, Hughes, 
Rakowski, & Decourchelle 2000). 



When the acceleration is efficient, the non-thermal 
population is expected to modify directly the ionization 
rates in the plasma as well as the line excitation (e.g. 
Dzifcakova |200C| , Seely et al. |1987| ). A hybrid electron 
distribution (Maxwellian plus power-law tail) is expected 
from diffusive shock acceleration (e.g., Berezhko & Ellison 



1999, Bykov & Uvarov 199E). The low energy end of the 



power-law electron distribution (which connects to the 
Maxwellian thermal population) is likely to enhance the 
ionization rates and to significantly modify the degree of 
ionization of the plasma, which is used as a diagnostic of 
the plasma electron temperature. 
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Fig. 1. The Hybrid (Maxwell/Power-law) electron distri- 
bution for different values of the break parameter x\> , com- 
pared to a Maxwellian distribution (xb = oo). The slope 
of the power-law has been fixed to a = 2. Reduced energy 
coordinates, x — E/kT, are used. 

In this paper, we shall examine the influence of a 
power-law non-thermal electron distribution (connecting 
to the falling Maxwellian thermal population) on the ion- 
ization and recombination rates for C, N, O, Ne, Mg, Si, 
S, Ar, Ca, Fe and Ni. For different characteristic values 
of the power-law electron distribution, the mean electric 
charge of these elements has been determined as a func- 
tion of the temperature at ionization equilibrium and for 
different values of the ionization timescale. We intend, in 
this paper, to give a comprehensive study of the depen- 
dence of these quantities on the parameters of the non- 
thermal population, illustrated by simple examples. We 
do not provide tables, which would be too numerous as 
the ionization equilibrium depends in our model on four 
parameters (element, temperature of the thermal compo- 
nent, index and low energy break of the non-thermal pop- 
ulation). In the appendix or directly in the text, we give 
the formula needed for the calculations of the rates which 
could be easily inserted in computer codes. 

In Section ^, we define the Hybrid electron distribution 
used in this work. The calculation of the new ionization 
collisional rates and (radiative and dielectronic) recombi- 
nation rates is discussed in Section ||. In Section [|, we 
present the derived mean electric charge of the elements 
in ionization equilibrium as well as in ionizing plasmas. 

2. The electron distribution shapes 

The Maxwellian distribution, generally considered for the 
electron distribution in astrophysical plasmas, N C (E), is 
defined as: 



dN e (E) = n c fg(E) dE 
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Fig. 2. Variation of the median energy (top panel) and 
of the normalisation constant (bottom panel) with the 
break parameter Xh- The lines are labeled by the value 
of the slope parameter a. The bold part of the curves 
corresponds to > a + 1/2. The black horizontal lines 
correspond to the values obtained for a Maxwellian distri- 
bution. 



where E is the energy of the electron, T is the electronic 
temperature and n e the total electronic density. In this 
expression the Maxwellian function f^(E) is normalised 
so that J °° /*?(£?) dE = 1. 

It is convenient to express this distribution in term of 
the reduced energy x = E/kT: 



dN c (x) 
A*) 



n c f M (x) dx 
2 



x 2 e 



(3) 
(4) 



The corresponding scaled (non-dimensional) distribution 
f M (x) is an universal function, of fixed shape. 

Non-Maxwellian electron distributions expected in the 
vicinity of shock waves, as in young supernova remnants, 
seem to be reasonably described by a Maxwellian distri- 
bution at low energy up to a break energy E^, and by a 
power-law distribution at higher energy (e.g., Berezhko & 
Ellison |1999| , Bykov & Uvarov |1999| ) . We call hereafter this 
"Maxwellian/Power-law" type of electron distribution the 
Hybrid electron distribution (/ H ). It is defined, in re- 
duced energy coordinates, as: 



dN e {x) = n c f n (x)dx 



(5) 
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f H (x) = C(x h ,a) — x 1 ' 2 e" 



x < Xb (6) 



/ H (x) = C(x h ,a) 4= x^ 2 f— ) z > x h , 

where Xb = E^/kT is the reduced break energy, and a 
is the energy index of the power-law (a > 1). Note that 
for a < 2 the energy diverges (in practice a cutoff at 
very high energy occurs). Since for the calculations of the 
ionization and recombination rates the very high energies 
(> 20kT) have negligible effect, for simplicity, we use here 
a power-law defined from Xb to infinity. 
The normalisation factor of the power-law distribution is 
defined so that the electron distribution is continous at x\, . 
The factor C(xb,ct) is a normalisation constant, so that 
V f K (x)dx 1: 



C(x h ,a) 



7 (§,x b ) + (a-l)- 1 



3/2 



(7) 



where 7(0, x) is the gamma function defined as 
7(0, x) = J Q t a ~ Y e _t dt. For x< Xb, the Hybrid distribu- 
tion only differs from a Maxwellian distribution by this 
multiplicative factor. The scaled distribution / H (x) only 
depends on the two non-dimensional parameters, Xb and 
a. The dependence on kT of the corresponding physical 
electron distribution is /£(£) = (kT)" 1 / H (£/kT). 

The Hybrid distributions / H (x), obtained for sev- 
eral values of the energy break Xb, are compared to the 
Maxwellian distribution in Fig. [j]. The slope has been 
fixed to a = 2, a typical value found in the models refer- 
enced above. The variation of the reduced median energy 
of the distribution with Xb-, for a — 1.5, 2., 3., is plotted in 
Fig. ||, as well as the variation of the normalisation factor 
C(x h ,a). 

As apparent in the figures, there is a critical value 
of Xb, for each a value, corresponding to a qualitative 
change in the behavior of the Hybrid distribution. This 
can be understood by looking at the distribution at the 
break energy Xb- Whereas the distribution is continuous, 
its slope changes. The logarithmic slope is 1/2 — Xb on the 
Maxwellian side and —a on the power-law side. There 
is no break in the shape of the Hybrid distribution (full 
line in Fig. [I]), only for the critical value of Xb = a + 1/2. 
For Xb > a + 1/2, the power-law always decreases less 
rapidly with energy than a Maxwellian distribution and 
does correspond to an enhanced high energy tail. The 
contribution of this tail increases with decreasing Xb 
(and a). Thus the median energy increases and the 
normalisation parameter, which scales the Maxwellian 
part, decreases (see Fig. |^). On the other hand, when 
Xb < a + 1/2, there is an intermediate region above 
the energy break where the power-law decreases more 
steeply than a Maxwellian (see dotted line in Fig. [j]). 
This results in a deficit of electrons at these energies 
as compared to a Maxwellian distribution, more and 
more pronounced as Xb is small. The median energy 



thus starts to decrease with decreasing Xb and can be 
even lower than the median energy of a Maxwellian 
distribution (Fig. ||). In this paper we only consider the 
regime where Xh > a + 1/2. It corresponds to clear cases 
where the high energy part of the distribution is indeed 
increased, as expected when electron are accelerated 
in shocks. Furthermore, the distribution used here is 
only an approximation, valid when the hard tail can be 
considered as a perturbation of the original Maxw ellian 
distribution. The simulations of Bykov and Uvarov ( |l999| , 
see their Figure 2) clearly show that the low energy part 
of the distribution is less and less well approximated by 
a Maxwellian distribution, as the 'enhanced' high energy 
tail extends to lower and lower energy (lower 'break'). 
Although we cannot rigorously define a corresponding 
quantitative lower limit on Xb, the cases presented by 
Bykov and Uvarov ( 1999 , see their Figure 2) suggest a 
limit similar to the one considered here, i.e. a few times 
the Maxwellian peak energy. 

The distribution considered here differs from the so- 
called "kappa-distribution" or the "power distribution", 
relevant for other physical conditions (see e.g. Dzifcakova 
2000 and references therein). These two distributions have 
been used to model deviations from a Maxwellian distri- 
bution caused by strong plasma inhomogeneities, as in the 
solar corona, and their impact on the ionization balance 
has been extensively studied (e.g. Roussel-Dupre 1980, 
Owocki & Scudder 1982, Dzifcakova 1992, Dzifcakova 



1998). Although the effect of the Hybrid distribution is ex- 
pected to be qualitatively similar, it has never been quan- 
titatively studied. In the next section we discuss how the 
ionization and recombination rates are modified, as com- 
pared to a pure Maxwellian distribution, depending on the 
parameters Xb and a. 

3. Calculations of the collisional ionization and 
recombination rates 

Let us consider a collisional process of cross section 
&{E), varying with energy E of the incident electron. 
The corresponding rate coefficient (cm 3 s _1 ), either for a 
Maxwellian distribution or a Hybrid distribution, f(x), is 
given by: 



Rate 



/2kT 
\ m 



x 2 cr(xkT) f(x) dx 



(8) 



with x t h = Eth/kT. E t b corresponds to the threshold en- 
ergy of the considered process (for E < E t b, cr(E) = 0). 
For the recombination processes, no threshold energy is 
involved and Xth = 0. 

The rates for the Hybrid distribution depend on kT, 
Xb and a and are noted C : H (T, Xb,ct), o^ R (T, Xb, a) and 
a r>R Xh ' a ) f° r ^ ne ionization, radiative and dielectronic 
recombination process respectively. The corresponding 
rates for the Maxwellian distribution which only depends 
on kT are C : M (T) , a^ R (T) and (T) . 
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Fig. 3. Variation of the ratio of the ionization rate in a 
Hybrid distribution over that in a Maxwellian with the 
same temperature, versus the break parameter of the 
Hybrid distribution. The curves correspond to different 
values of the slope parameter a and are labeled accord- 
ingly. The ion considered is +6 , the temperature is fixed 
at 10 6 K. 



The ionization data are taken from Arnaud & 
Rothenflug ( |1985| ) and Arn aud fc Raymond ( |1992| ), as 
adopted by Mazzotta et al. ( |1998| ) for the most abundant 
elements considered here. The recombination data are 
taken from the updated calculations of Mazzotta ( 1998| ) . 
In the next sections we outline the general behavior of 
the rates with the electron distribution parameters, using 
mostly oxygen ions (but also iron) as illustration. 



3.1. The electronic collisional ionization rates 

The ionization cross sections present a threshold at the 
first ionization potential of the ionizing ion, E I . The cross 
sections always present a maximum, at E m , and decrease 
as \n(E)/E at very high energies (e.g., Tawara, Kato,& 
Ohnishi 1985). The ionization rate is very sensitive to the 
proportion of electrons above the threshold and the mod- 
ification of the ionization rate for the Hybrid distribution 
depends on how the high energy tail affects this propor- 
tion. 

Parametric formulae for the ionization cross sections 
are available from the litterature and it is easy to derive 
the corresponding rates for the Hybrid distribution. This 
is detailed in Appendix [A|. 

To understand the influence of the presence of a high 
energy power-law tail in the electron distribution, we 
computed the ratio ^(T, x h , a) = C : H (T, x^, a)/C : M (T), 
of the ionization rate in a Hybrid distribution over that 
in a Maxwellian with the same temperature. This ratio is 
plotted in Fig. || to Fig. |?j for different ions and values of 
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Fig. 4. Same as Fig. ^ but the parameter a is now fixed 
to 2. The curves correspond to different values of the tem- 
perature. 
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Fig. 5. Same as Fig. |3| for the different ions of oxygen. 
Each curve is labeled by the charge of the ion considered. 
The slope parameter is fixed to a — 2 and the temperature 
for each ion is fixed at the value, T* , where the abundance 
of the ion is maximum for a Maxwellian electron distribu- 
tion. 



the parameters Xh and a. 

Let us first consider +6 . Its ionization potential 
is Ei — 739 eV and the cross section is maximum at 
about 3 E\. Its abundance, for a Maxwellian electron 
distribution, is maximum at T* ~ 10 6 K under ionization 
equilibrium (Arnaud & Rothenflug 1985). At this tem- 
perature, the threshold energy is well above the thermal 
energy (E x /kT ~ 8) and only the very high energy tail of 
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Fig. 6. Same as Fig. [| for the different ions of iron. 

the Maxwellian contributes to Cj M (T), i.e. a small fraction 
of the electron distribution. This fraction is dramatically 
increased in the Hybrid distribution as soon at the break 
energy is not too far off from the threshold, Xb ~ 15 
for +6 (Fig. ||). The enhancement factor /3 : (T, Xb,a) 
naturally increases with decreasing break a; b and slope a 
parameters (Fig. since the distribution median energy 
increases when these parameters are decreased (Fig. ^|). 

This behavior versus Xb and a is general at all 
temperatures as illustrated in Fig. [|, provided that the 
thermal energy is not too close to E m , i.e. that the 
majority of the contribution to the ionization rate is from 
electrons with energies corresponding to the increasing 
part of the ionization cross section. If this is no more the 
case, the ionization rate starts to decrease with increasing 
distribution median energy. Thus, for high enough values 
of the temperature (see the curve at T = 10 s K in 
Fig. |J), the factor /^(T, x\,,a) becomes less than unity 
and decreases with decreasing Xb- The correction factor 
is small (around ~ 10%) however in that case. 

More generally the enhancement factor P T (T, Xh,ct) 
at fixed values of Xb and a, depends on the temperature 
(Fig. |^). It decreases with increasing temperature: the 
peak of the distribution is shifted to higher energy as the 
ratio kT/E l increases and the enhancement due to the 
contribution of the hard energy tail decreases. 

The qualitative behavior outlined above does not 
depend on the ion considered. We plotted in Fig. [| 
and in Fig. || the enhancement factor for the different 
ions of oxygen and a choice of iron ions at T* (the 
temperature of maximum ionization fraction of the ion 
for a Maxwellian electron distribution under ionization 
equilibrium) . E Y /kT* is always greater than unity and the 
ionization rates are increased by the Hybrid distribution, 



the enhancement factor /3 t (T, x\>,a) increasing with 
decreasing Xb- However this enhancement factor differs 
from ion to ion, it generally increases with increasing 
E l /YT* value (approximatively with an exponential 
dependence), as shown in Fig. [?]. This is again due to 
the relative position of the peak of the distribution with 
respect to the threshold energy. Note that Ej/kT* is 
generally smaller for more ionized ions (but this is not 
strictly true) so that low charge species are generally 
more affected by the Hybrid distribution. 

In summary, the Hybrid rates are increased with 
respect to the Maxwellian rates except at very high 
temperature. The enhancement factor depends on the 
temperature, mostly via the factor E l /YT. It increases 
dramatically with decreasing temperature and is always 
important at T*, where it can reach several orders of 
magnitude. The ionization balance is thus likely to 
be affected significantly, whereas the effect should be 
smaller in ionizing plasmas but important in recombining 
plasmas. For Xb typically lower than 10-20 (with this 
upper limit higher for lower temperature, see fig. ^), the 
impact of the Hybrid rate increases with decreasing Xb 
and a. 

The ionization rates for a Hybrid distribution are less 
dependent on the temperature than the Maxwellian rates, 
as illustrated on Fig. | and on Fig. § This is a direct 
consequence of the temperature dependance of the en- 
hancement factor: as this factor increases with decreas- 
ing temperature, the Hybrid ionization rate decreases less 
steeply with temperature than the Maxwellian rates. More 
precisely, as derived from the respective expression of the 



rates at low temperature (respectively Equation A. 2 and 
Equation A. 10), the Maxwellian rate falls off exponen- 
tially (as c~ K i /kT ) with decreasing temperature, whereas 
the Hybrid rate only decreases as a power-law. As ex- 
pected, one also notes that the modification of the rates is 
more pronounced for lower value of Xb (compare the two 
figures corresponding to Xb — 10 and Xb = 2.5). 

3.2. The recombination rates 

3.2.1. The radiative recombination rates 

The radiative recombination rates are expected to be 
less affected by the Hybrid distribution, since the cross 
sections for recombination decrease with energy and no 
threshold exists. As the net effect of the high energy 
tail present in the hybrid distribution is to increase the 
median energy of the distribution (cf. fig. ||) , as compared 
to a Maxwellian, the radiative recombination rates are 
decreased. 

To estimate the corresponding dumping factor, 
(3 RR (xb, a) — a RR (T, Xb, a)/ a^ ( T), we follow the method 
used by Owocki & Scudder (1983). We assume that the ra- 
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diative recombination cross section varies as a power-law 
in energy: 

a RR (E)^E- a (9) 

which corresponds to a recombination rate (Equation ^J), 
for a Maxwellian distribution (Equation |J), varying as: 

a M (T) cx T v 

with i] = a — h . 



(10) 



The dumping factor computed for such a power-law 
cross section (Equation || with Equation ^|) is: 



Am Ob, a) 



f^x-v / H (x) dx 
f^x-v f M (x) dx 



(11) 



Note that the dumping factor is independent of the 
temperature. It depends on the ion considered via the 
rj parameter. Replacing the Maxwellian and Hybrid 
distribution functions by their expression (respectively 
Equations ^ and ||) we obtain: 



Am Ob,") 



C(x h ,a) 



^2 



r),x h 



(12) 



This estimate of the dumping factor is only an approx- 
imation, since the radiative recombination has to be com- 
puted by summing over the various possible states of the 
recombined ions, taking into account the respective differ- 
ent cross sections. Furthermore, even if often the radiative 
recombination rate can be approximated by a power-law 
in a given temperature range, this does not mean that the 
underlying cross section is well approximated by a unique 
power-law. 

However as we will see the correction factor is small, 
and we can reasonably assume that it allows a fair esti- 
mate of the true Hybrid radiative recombination rates. To 
minimize the errors, the Hybrid radiative recombination 
rate has to be calculated from the best estimate of the 
Maxwellian rates, multiplied by this approximation of the 
dumping factor: 



P KK (x h ,a) cgjT) 



(13) 



where a^f B (T) is as given in Mazzotta et al. (|199SD . 

The parameters r\ for the v arious ions are taken from 
Aldrovandi & Pequignot ( 1973 ), when available. For other 
ions we used a mean value of r\ = 0.8 corresponding to 
the mean value < ?y > reported in Arnaud & Rothcnflug 
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( |l985|) . The exact value has a negligible effect on the 
estimation of the radiative recombination rates. 
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Fig. 10. Dumping factor for the radiative recombination 
rates, (3 KR (xb,oe), of oxygen ions versus the break pa- 
rameter Xb- The curves correspond to different values of 
the slope parameter a and are labeled accordingly. This 
dumping factor is independent of temperature. 
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The slight variations among the species are due to the 
corresponding variations in the adopted 7/ values (see 
Equation O) . 



The dumping factor is plotted in Fig. |l0| for the various 
ions of oxygen. In that case a common rj value is used. The 
dumping factor decreases with decreasing values of x\> and 
a, following the increase of the distribution median energy. 
The modification is however always modest, at most 15% 
for a = 2. For iron, plotted in Fig. [ll] for a — 2, the value 
of r\ slightly changes with the considered ions, but this 
only yields negligible variations in the dumping factor. 



3.2.2. The dielectronic recombination rates 

The dielectronic recombination is a resonant process in- 
volving bound states at discrete energies Ei and the 
rates have to be computed by summing the contribu- 
tion of ma ny su ch bound states. Accor ding t o Arnaud & 
Raymond ( 1992| ), and Mazzotta et al. ( 1998 ), the dielec- 
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Fig. 12. Variation of the ratio of the DR rate in a Hybrid 
distribution over that in a Maxwellian with the same tem- 
perature, versus the break parameter Xb of the Hybrid 
distribution. The curves correspond to different values of 
the slope parameter a and are labeled accordingly. The 
ion considered is +6 , the temperature is fixed at 10 6 K. 

tronic recombination rates for a Maxwellian distribution 
can be fitted accurately by the formula: 



, M 



(T)=T - v 3/2 e _„ 



3 _ i 

cm s 



(14) 



where T c y is the temperature expressed in eV and 
X{ — Ei/kT. The numerical values for a and E\ are 



taken from Mazzotta et al. (1998). Only a few terms 
(typically 1 to 4) are introduced in this fitting formula. 
They roughly correspond to the dominant transitions for 
the temperature range considered. 

Fo llowing again the method used by Owocki & Scudder 
(|1983|) , we thus assume that the corresponding dielec- 
tronic recombination cross section can be approximated 
by: 



d 5(E-Ei) with Ci = 



(27rm c )' 
4 Ek 



(15) 



The relation between C\ and c; is obtained by comparing 
Equation (nj) with the equation obtained by integrating 
(Equation g) the above cross section over a Maxwellian 
distribution (Equation ||) . The dielectronic rates can then 
be computed from Equation (||), with the cross section 
given by Equation ( |l5| ) and the distribution function given 
by Equation (^|): 



o^ R (T,x h ,a) = C(x h ,a) T cV 3/2 ^ a e Xi 
+ C(x h ,a) r c " 3/2 e~^ ]T a 



(16) 
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Fig. 13. Same as Fig. |l2] but the parameter a is now fixed 
to 2. The curves correspond to different values of the tem- 
perature and are labelled accordingly. 

Note that this estimate of a^ R (T, xt,, ot) is only an 
approximation, for the same reasons outlined above for 
the radiative recombination rates. 

To understand the effect of the hybrid distribution, 
let us assume that only a single energy E DR is dominant, 
corresponding to a simple Dirac cross section at this en- 
ergy. In that case, from Equation (ft), the ratio of the 
dielectronic recombination rate in a Hybrid distribution 
over that in a Maxwellian with the same temperature, 
(3 nli {T,x h ,a) = a^ R (T,x h ,a)/a^(T), is simply the ra- 
tio of the Hybrid to the Maxwellian function at the res- 
onance energy. Its expression depends on the position of 
the resonance energy with respect to the energy break. In 
reduced energy coordinates, we obtain from Equation (jij) 
and Equation (0): 



C(x h ,a) 



C(x\>, a) e 



kT 



< x h 



(17) 



kT 



XbkT 
> x h 



-(!+«) 



For (xb kT) > E Dn (i.e. at high temperature or high 
value of Xb), the resonance lies in the Maxwellian part of 
the distribution. f3 DR (T, Xb, a) is independent of the tem- 
perature and the dielectronic recombination rates are de- 
creased, following the variation of the normalisation fac- 
tor, C(xb,a), i.e. the decrease is modest (see fig. ||). 

For (xb kT) < E DR the resonance lies in the power-law 
part of the distribution. The increase of the dielectronic 
recombination rate can be dramatic, increasing with 
decreasing Xb and a. 
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These effects of the Hybrid distribution on the dielec- 
tronic recombination rates are illustrated in Fig. |l2| to 
Fig. [l5|, where we plotted the factor /3 DR (T, Xb, a) for var- 
ious ions and values of the parameters. The factors are 
computed exactly from Equations (Q) and (jl6|). 

In Fig. [l^ we consider +6 at the temperature of 
its maximum ionization fraction, T* = 10 6 K. For this 
ion only one term is included in the rate estimate, with 
E DT< = 529 eV, and £' DR /kT = 6.1 at the temperature 
considered. We plotted the variation of [3 DR (T, Xb, a) with 
£b for a — 3, a = 2 and a = 1.5. For Xb > 6.1 the 
'resonance' energy E OR lies in the Maxwellian part of the 



distribution and the dielectronic recombination rate is de- 
creased as compared to a Maxwellian, but by less than 
10%, following the variation of the normalisation factor 
C(xb,a). For smaller values of Xb, the rate is increased 
significantly, up to a factor of 5 for a = 1.5. 

We consider other temperatures, fixing a to a = 2, 
in Fig. 13. Since we only consider the parameter range 
x\, > a + 1/2, there is a threshold temperature, kT > 
E UR / (a + 1/2), above which the resonance always falls in 
the Maxwellian part. The dielectronic recombination rate 
is decreased via the factor C(xb,a). This factor slightly 
decreases with decreasing Xb (cf . fig. ||) . At lower temper- 
ature, the resonance energy can fall above the break, pro- 
vided that Xb is small enough (xb < E DR /YT). This occurs 
at smaller Xb for higher temperature and the enhancement 
at a given Xb increases with decreasing temperature. 

We display in Fig. [l4| and Fig. [l5] the variation of 
the factor /3 DR (T, Xb,ot) with Xb (for a = 2), for the 
different ions of oxygen and iron, at the temperature 
of maximum ionization fraction for a Maxwellian dis- 
tribution under ionization equilibrium. For most of the 
ions this temperature is above the threshold temperature, 
kT = E- DIi /(a + 1/2), for all the resonances and the di- 
electronic rate is decreased. For the ions for which this 
is not the case (0 +1 , 0+ 6 and from Fe +1 to Fe +5 ), the 
dielectronic rate can be increased significantly (by a fac- 
tor between 2 to 5) provided Xb is small enough (typically 
Xb = 2.5—5). The increase starts as soon as Xb < E OR /YT* 
for the oxygen ions. The behavior of /3 DR (T, Xb, a) is more 
complex for the iron ions (two breaks in the variation of 
/3 DR (T, Xb, a)), due to the presence of more than one dom- 
inant resonance energy (more than one term) , taken into 
account in the computation of the dielectronic rate. 

In conclusion, the effect of the hybrid distribution on 
the dielectronic rate depends on the position of the reso- 
nance energy as compared to the power-law energy break. 
It can only be increased if kT < E DB /(a + 1/2). At 
high temperature, the dielectronic recombination rate is 
slightly decreased. 

3.2.3. The total recombination rates 

At Xb = 10, the total rates are basically unchanged by the 
Hybrid distribution. For x h = 2.5 = a+ 1/2 (Fig. |l|j, the 
total rates are more significantly changed. The radiative 
recombination rate increases with decreasing temperature 
and it usually dominates the total recombination rate in 
the low temperature range. As the dielectronic rate is in- 
creased by the Hybrid distribution only at low tempera- 
ture, there are very few ions for which the total recombi- 
nation rate can be actually increased. This only occurs in 
a small temperature range, in the rising part of the dielec- 
tronic rate. One also notes the expected slight decrease 
of the radiative recombination rates (when it is dominant 
at low temperature) and of the dielectronic rate at high 
temperature. 



10 



Porquet, Arnaud, Decourchelle: Impacts of a power-law electron tail 




Temperature [Kelvin] 



Fig. 16. Variation of the total recombination rate with 
temperature. Each curve corresponds to an ion and is la- 
beled accordingly. Black lines: rates for a Hybrid electron 
distribution with a = 2 and Xb = a+ i = 2.5. Black thick 
lines: rates for Maxwellian distribution. 

4. Ionization equilibria 

4.1. Collisional ionization equilibrium (CIE) 

The ionization equilibrium fractions, for coronal plasmas, 
can be computed from the rates described in the previous 
sections. In the low density regime (coronal plasmas) the 
steady state ionic fractions do not depend on the electron 
density and the population density ratio Nz.z+i/Nz, z of 
two adjacent ionization stages Z + ( z+1 ^ and Z +z of element 
Z can be expressed by: 

N z,z+i _ Cf' z 

~^z7 - (18) 

where Cj Z,z and o^ z+1 are the ionization and total 
recombination rates of ion Z +z and Z + ^ z+1 ^ respectively. 
To assess the impact of the Hybrid rates on the ionization 
balance, we computed the variation with temperature of 
the mean electric charge of the plasma. This variation is 
compared with the variation obtained for a Maxwellian 
electron distribution on Fig. ^ for oxygen and iron and 
for different values of the parameters Xb and a. 

As expected, the plasma is always more ionized for a 
Hybrid electron distribution than for a Maxwellian dis- 
tribution. The mean charge at a given temperature is in- 
creased, since the enhancement of the ionization rate is 
always much more important than a potential increase of 
the dielectronic rate (e.g. compare Fig. | and Fig.|l|). The 
effect of the Hybrid distribution on the plasma ionization 
state is thus governed by the enhancement of the ioniza- 
tion rates. The enhancement of the plasma mean charge is 
more pronounced for smaller values of Xb and smaller val- 



ues of a (Fig. [n]) , following the same behavior observed for 
the ionization rates (due to the increasing influence of the 
high energy tail). Similarly the effect is more important 
at low temperature, and a clear signature of the Hybrid 
distribution is the disappearance of the lowest ionization 
stages, that cannot survive even at very low temperature. 
For instance, for a = 2 and the extreme corresponding 
value of Xb = a + 1/2, the mean charge is already +4 
for oxygen and +6 for iron at T — 10 4 K. At high tem- 
perature, the mean charge can typically be changed by a 
few units, the effect being more important in the temper- 
ature range where the mean charge changes rapidly with 
temperature in the Maxwellian case. 

The same behavior is seen for all elements (Fig. |l8| ). 
One notes that the effect of the Hybrid distribution 
generally decreases with Z. Again this is a consequence of 
the same behavior observed on the ionization rates (see 
Fig.0). 

A remarkable effect of the Hybrid distribution is that 
the mean charge is not always a monotonous function 
of temperature, in the low temperature regime. This is 
clearly apparent in Fig. || and Fig. [l| for 10 4 K < T < 
10 5 K and Xb = 10. This phenomenon can only occur 
when the dielectronic rate dominates the total recombi- 
nation rate and in the temperature range where this rate 
increases with temperature. In that case, the density ratio 
of two adjacent ions, A^z, z +i/ Nz,z, can decrease with tem- 
perature provided that the ionization rate of Z +z increases 
less rapidly with temperature than the recombination rate 
of the adjacent ion (EqJI"g|). This usually does not 

occur in the Maxwellian case, but can occur in the Hybrid 
case, due to the flatter temperature dependence of the ion- 
ization rates for this type of distribution. For instance, for 
3 10 4 K < T < 7 10 4 K, the ionization rate of 0+ 2 is in- 
creased by a factor of 2.5 for an Hybrid distribution with 
Xb = 10 (Fig. ||), whereas the total recombination rate 
of +3 is increased by a slightly larger factor of 2.7 (see 
the corresponding grey line in Fig. as seen above for 
Xb=10 the total rate is basically unchanged compared to 
the Maxwellian case). The mean charge, which is around 
(z) — 2.5, is thus dominated by the behavior of these ions 
and decreases in that temperature range. 

4.2. Non-equilibrium ionization (NEI) 

Collisional Ionization Equilibrium (CIE) is not always 
achieved. For example, in adiabatic supernova remnants, 
the ionization timescale is longer than the dynamical 
timescale, so that the plasma is underionized compared to 
the equilibrium case. In non-equilibrium conditions, the 
ionization state of the gas depends on the thermodynamic 
history of the shocked gas (temperature, density) and time 
elapsed since it has been shocked. 
The time evolution of the ionic fractions is given by: 

u = ricl^' A ZiZ _i + a R ' A ZjZ+ i (19) 
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Fig. 17. Mean electric charge versus temperature for oxygen (top panels) and iron (bottom panels) for different electron 
distributions. Black thick lines: Maxwellian electron distribution. Black thin lines: Hybrid distribution with a = 1.5 
(left panels), a — 2 (middle panels) and a — 3 (right panels). Three values of Xb are considered: = 10 (dotted line), 
Xb = 5 (dashed line) and the extreme value Xb = a + 1/2 (full line). 



a^ z )X z , z ] 



with Xz z = 



N, 



Z.z 



Z.i 



To estimate the effects of a Hybrid electron distribution 
on the ionization in non-equilibrium ionization conditions, 
we assume that the gas has been suddently heated to a 
given temperature, which stays constant during the evo- 
lution. The ionization timescale depends then on J n e dt, 
where n c is the number density of electrons and t the 
time elapsed since the gas has been heated. Within this 
assumption, the coupled system of rate equations can be 
resolved using an exponentiation method (e.g., Hughes & 
Helfand 1985). 

For different ionization timescales (up to equilibrium) , 
we computed the variation with temperature of the mean 
electric charge of oxygen and iron in two extreme cases 
of the electron distribution: Maxwellian and Hybrid with 
Xb = a + 1/2 and a = 2. For small ionization timescales 
(n e t ~ 10 s — 10 9 s cm -3 ), the effect of the Hybrid distribu- 



tion on the mean electric charge is small, it increases with 
the ionization timescale and is maximum at equilibrium 
as is illustrated for oxygen and iron in Fig. 19. As in the 



equilibrium case, the effect from non-thermal electrons is 
always more important at low temperature and vanishes 
at high temperature. Note that the mean electric charge 
is slightly larger at high temperature for the thermal pop- 
ulation than for the non-thermal one, as a consequence of 
the decrease of the ionization cross section at very high 
energy. 

5. Conclusions 

We have studied the effect on the ionization and recom- 
bination rates, as well as on the ionization balance, of 
a non-thermal electron distribution, as expected in the 
vicinity of strong shocks. 

The electron distribution is modelled by a Maxwellian 
distribution at low energy up to a break energy, and by 
a power-law distribution at higher energy. It is carac- 
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Fig. 18. Mean electric charge versus temperature for the Hybrid electron distribution (black thin lines) compared to 
the Maxwellian distribution (black thick lines). The slope parameter a of the Hybrid distribution is fixed to 2, and 
the Xb values are as in Fig. |l^. Each panel corresponds to an element and is labelled accordingly. 



terised by the three parameters kT (the temperature of 
the Maxwellian part), Xh the reduced energy break, and 
a the slope of the power-law component. We only consid- 
ered the parameter range where Xb > a + 1/2 which corre- 
sponds to an enhanced high energy tail. All the behaviors 
outlined are only valid for this range of parameters. 

We provide exact formulae of the ionization rates for 
this Hybrid electron distribution in the Appendix, and ap- 
proximate estimates of the radiative recombination rates 
(Eq. [12] and Eq. |l^) and of the dielectronic recombination 
rates (Eq. |l^). The Hybrid rates depend on the ion con- 



sidered and on the parameters kT, a and Xh- Computer 
codes are available on request. 

For the parameter range considered, the proportion 
of electrons at high energies and the mean energy of the 
distribution is a monotonic function of Xb and a. As ex- 
pected, the modification of the rates for the Hybrid dis- 
tribution, as compared to the Maxwellian distribution of 
the same temperature, increases with decreasing x\, (with 
a threshold at about x\, ~ 10 — 20, higher for lower tem- 
perature) and decreasing a. 

The impact of the Hybrid electron distribution on the 
ionization rates depends on how the high energy tail af- 
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Fig. 19. Mean electric charge versus temperature for oxygen (top panels) and iron (bottom panels) for different 
ionization timescales (thin lines) compared to equilibrium (bold lines) and for two extreme electron distributions, 
namely Maxwellian (solid lines) and Hybrid (dotted lines) for a = 2 and Xb = a + 1/2. 



fects the proportion of electrons above the ionization po- to the Maxwellian rates except at very high tcmpera- 
tential E I . The Hybrid rates are increased with respect ture. The enhancement factor depends on the tempera- 
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ture, mostly via the factor E 1 /YT, and increases dramat- 
ically with decreasing temperature. For a given ion, it is 
always important at T*, the temperature of maximum ion- 
ization fraction for a Maxwellian distribution under ion- 
ization equilibrium, where it can reach several orders of 
magnitude. 

The effect of the hybrid distribution on the diclectronic 
rate depends on the position of the resonance energies E DR 
as compared to the power-law energy break. The dielec- 
tronic rate can only be increased if kT < E DR /(a + 1/2). 
At T* the enhancement factor is typically less than an 
order of magnitude. At high temperature, the dielectronic 
recombination rate is slightly decreased (by typically 10% 
at most). The effect of the hybrid distribution on the ra- 
diative recombination rates is only of the order of a few 
10% at most. 

The ionization balance is affected significantly, whereas 
the effect is smaller in ionizing NIE plasmas. The plasma 
is always more ionized for a Hybrid electron distribution 
than for a Maxwellian distribution. The effect is more im- 
portant at low temperature, and a clear signature of the 
Hybrid distribution is the disappearance of the lowest ion- 
ization stages, which cannot survive even at very low tem- 
perature. 
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Appendix A: Ionization rates 

A.l. Direct ionization (Dl) 

For the direct ionization (DI) cross sections we chose the 
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For Xb > xy. 

(xj, Xh, en) is the sum of the contribution of the 
truncated Maxwellian component and the power-law 
component: 



fitting formula proposed by Arn aud & Rothenflug (1985) 
from the work of Younger ( 1981 ): 



G*(x h x h ,o) 



-X 



F {x\ ) 

DI v j / 



1 J .1 



Ai C/j + Bj C/f + Cj Info) ■ lh 



ln(i 



x i x i 



where: 



with u\ = — 



U, = l- 



(A.l) 

The sum is performed over the subshclls j of the ionizing 
ion. E is the incident electron energy and 7j is the colli- 
sional ionization potential for the level j considered. 
The parameters Aj, Bj, Cj, Dj (in units of 10~ 14 cm 2 eV 2 ) 
and Ij (in eV) are taken from the works of Arnaud & 
Raymond (1992) for iron, and of Arnaud & Rothenflug 
(1985) for the others elements. The parameters for ele- 
ments not considered in these works are given in Mazzotta 
et al. ( p98| ). 

A. 1.1. The Maxwellian electron distribution 

For a Maxwellian electron distribution, Arnaud & 
Rothenfl ug (|1985 ) obtained according to Equations (|l|), 
(§) and JAT), the rate: 



Xh 

Wbj = — 



(A.6) 
(A.7) 



^ni^j^b) = Aj [l-Xj/l(x b )] 



l + -i--a; i (2 + a: j )/i(xb) 



Cj [/i(x b )+ln( Wbj )] (A.S 
Dj [xjf 2 (xb) + hi (ttbj) Xj/i(x b )] 

"bj 1 

a -1/2 a + 1/2 

"bj _ 2 < 
a -1/2 a + 1/2 a + 3/2 



.4 



Si 



C M (T) 

where 



6.692 x 10 7 

(kT) 3 / 2 



J2 e — F?M) cmV1 (A.2) 



a - 1/2 
a + 1/2 



ln(w bj ) 
ln(u bj } 



1 



a - 1/2 
1 



(A.9) 



a +1/2 



For x b < Xj: 

Only the power-law component contributes of the elec- 
tron distribution to the rate: 



1 3 

kT 



F DI (xj) = [l-xj/^xj)] 

+ Bj [l + Xj-Xj(2 + x j )/ 1 (x j )] 
+ Cj /i(xj) + TJj Xj/ 2 (xj) 



(A.3) 
(A.4) 



G^(xj,x b ,a) = 
where 

/di(") = 



-) / M («) 



2a 



(A.10) 



1/4 (a 2 -1/4) (a + 3/2) 



(a -1/2)' (a + 1/2)' 



where kT and 7j are in eV. The summation is per- 
formed over the subshells j of the ionizing ion. The math- 
ematical functions, /i(x) = e x ^-j— dt, and / 2 (x) = 

e * Ji° m W ^ can ^ e com P u ted from the a nalyt ical A. 2. Excitation autoionization (EA) 
approximations given by Arnaud & Rothenflug (1985) in 
their Appendix B. 



(A.ll) 



A. 1.2. The Hybrid electron distribution 

Similar to the Hybrid electron distribution, the direct ion- 
ization rate C5 :(T, Xb, a) is given by: 



For the excitation autoionization (EA) cross sections, 
we used the generalized formula proposed by Arnaud & 
Raymond ( p.992| ): 

v EA {E) = —— A + BU + CU 2 + DU 3 + E ln(u)] 



Ulr, 



with u 



E 



1 



1 



C D (T, ib, a) = C(x b ,a) 



6.692 x 10 7 

(kT) 3 / 2 



(A.12) 



^G"(xj,x b ,a) cm 3 s 1 



(A.5) 



where 7 EA is the excitation autoionization threshold and 
E is the incident electron energy. 

The parameters A, B, C, D, E (in units of KT 16 cm 2 eV) 
and I EA (in eV) a re taken from the works of Arna ud & 
The function (xj, x b , a) depends on the position of the Rothenflug ( 1985 ) and Arnaud & Raymond ( 1992 ). The 
power-law break energy as compared to the ionization po- parameters for elements not co nsidered in these works arc 
tential: given in Mazzotta et al. ( 1998 ). 
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A. 2.1. The Maxwellian electron distribution 

The excitation autoionization rate for a Maxwellian dis- 
tribution is: 



where 



6.692 x 10 7 e _I EA 

(kT)V2 



J EA 

kT 



+ C[l - x EA + xij^x^)} 



+ D 



1 - 



EA 
• l -E 

2 



(A.14) 
(A.15) 



+ EMx, 



A. 2. 2. The Hybrid electron distribution 

For the Hybrid electron distribution, the excitation au- 
toionization rate C E (T, x b , a) is given by: 



C*(T,x h ,a) = C(x h ,a) 



6.692 x 10 7 



+ D 
+ E 



1 



,-3 



a - 1/2 
1 



a + 5/2 
ln(u E 



+ 



(a -1/2) 2 a -1/2 



For x b < x EA : 

Only the power-law component contributes of the elec- 
tron distribution to the rate: 



G EA (x EA ,x h ,a) = x h e Xb [ 



(A.21) 



where 

/ea(") 



B 



a -1/2 



1/4 



+ 
+ 
+ 



2 C 



(a-l/2)(a + 3/2) 
3 £> 

(a-l/2)(a + 5/2) 



(A.22) 



(a -1/2)= 



cm 3 s 1 



A3. Tota/ ionization rates (Dl + EA) 

(A. 16) The total ionization rate C : H (T, Xb, a) is obtained by: 

(A.23) 



(kT)V2 
x G^ A (x EA ,a; b ,a) 

The function G EA (x EA , x b , a) depends on the position of Cf(T,Xh,a) = (T, Xb, a) + C EA (T, Xb, a 
the power-law break energy as compared to the ionization 
potential: 



- For x b > x EA : 

G E i A (x EA ,Xb,a) is the sum of the contribution of the 
truncated Maxwellian component and the power-law 
component: 



G* A (x EA ,x h ,a) = e-**A F™(x EA ) 

- ^ ^EA^EA^b) 



where: 



F EA (x EA ,x h ) = A + 5[l-x EA /i(x b )] 



+ C 



+ D 



l - 



"EA 

X 



(A.17) 



(A.18) 



(A.19) 



I v _|_ EA 



+ B 
+ C 



2«L ' 2m ba 
+ is[ln( UEA ) + / 1 (xb)] 
1 



- ~f /l(Xb) 



.a- 1/2 

1 "ea 1 

a -1/2 a + 1/2 



a -1/2 a + 3/2 



(A.20) 



